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Abstract. We are concerned with the first hitting times of the Bessel 
processes. We give explicit expressions for the densities by means of 
the zeros of the Bessel functions and show their asymptotic behavior. 
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1. Introduction 

The Bessel process is one of the fundamental stochastic processes. If the index is 
a half integer greater than or equal to —1/2, it is identical in law with the radial 
motion of a Euclidean Brownian motion. Moreover, in the black-Scholes model, 
a basic model in mathematical finance, a stock price processes is modelled by a 
geometric Brownian motion, which is a time change of some Bessel process. It is 
also known that the explicit computations on the Bessel process are useful in the 
study of so-called the CIR model of interest rates. 

The first hitting time to a point or a domain of stochastic process is an im- 
portant object in probability theory. When a geometric Brownian motion or a 
Bessel process is used to modell a stock price, the first hitting time plays an im- 
portant role in the theory of exotic options. For applications of Bessel processes 
to mathematical finance, see e.g., Geman-Yor [3] and Yor |12j . 

In this article we are concerened with the first hitting time to a point of the 
Bessel process. The purpose is to show an explicit expression for the density 
and, by using it, to show asymptotic behavior of the density at infinity. The 
probability distributon has been recently studied in [I] and we use the results 
therein. 

Let t^I be the first hitting time to b of a Bessel process with index v start- 
ing from a. When a < b, general theory on the eigenvalue expansions for the 
first hitting times of diffusion processes (cf. Kent [7]) may be applied since the 
boundary is not natural, and an explicit simple expression via the zeros of the 
Bessel function for the density of is known. However, in the case of a > b, 
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we cannot apply the general theory since the boundary oo is natural. Hence we 
concentrate on the latter case. 

By the theory of one-dimensional diffusion processes, the Laplace transform 
of the density (the moment generating function) of rf$ can be computed and is 
expressed as a ratio of the modified Bessel functions. We will invert the Laplace 
transform along the same line as that in |3], where we concentrate on the dis- 
tribution function. In this sence this article may be regarded as a companion of 

Recently Byczkowski et al [TJ [2] have obtained other expression for the density 
and have shown the order of its decay at infinity. In some cases (see Remark 13. 2p 
we need to use their results to show the exact asymptotics of the decay or to 
determine the constant. It should be also mentioned that Uchiyama [10J has 
studied the asymptotic behavior of the density in the case when v — 0. 



2. Density of the first hitting time 

For v G K the diffusion process on [0, oo) with infinitesimal generator 

^ M= l^_ + 2i/ + l d 



2 dx 2 2x dx 

is called the Bessel process with index v. If 2v + 2 is a positive integer, the 
Bessel process is identical in law with the radial motion of a (2i/ + 2)-dimensional 
Brownian motion. Hence the number 2v + 2 is called the dimension of the Bessel 
process. For all v the boundary oo is natural. The classification of depends on 
v, but it is not natural for any v. For details, see [SJ [6J, [T] . 

For a, b > 0, let denote the first hitting time to b of a Bessel process 
with index v starting from a. Then it is a fundamental fact in the theory of 
one- dimensional diffusion processes (cf. [5]) that, for A > 0, the function 

x i — y E[e *- b ] = u(x, A) 
is increasing (decreasing) on [0,6) (resp. (6, oo)) and satisfies u(b,\) = 1 and 
(2.1) ^ v) u = \u, 

where E denotes the expectation. Hence, solving this equation, we have an 
explicit expression for the Laplace transorm of the distribution of by means 
of the modified Bessel function. 
If < a < b, we have (see [6]) 

. At M. a- v I v (aV2X) 



Me T "' 1 ] = " v ' , v > -1, 

b-"I v (by/2X) 



and 



E[e-^] = a ~ VI -» {ay/ ™\ u<l, 
b-"I_ u {bV2X) 
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where l u is the modofied Bessel function. Since is not a natural boundary, the 
solution u(x, A) for the equation (12. ip may be written as the canonical product 
(see [Hi) 



u(x, A) = u(x 



on( 1+ £) 



where {A^}^ is a sequence of the simple and positive zeros of u(x, •). Therefore 
the distribution of may be written as an infinite convolution of mixtures of 
exponential distribution. For the general theory, see Kent [7]. 
On the other hand, if < b < a, we have 

2.2 E\e~ a > b = vK ' , i/ > -1, 

V ; L J b-"K v (bV2\) 

for every i/el, where K v is the other modified Bessel function called the Mac- 
donald function. The boundary oo is natural and the general theory above are 
not applicable. 

In the following we are only concerned with the latter case. We invert the 
Laplace transform (12. 2 p and we give an expression for the density of rf^ from 
Theorem 12.11 below, which is Theorem 3.3 in [I], on some ratio of the modified 
Bessel functions. 

We denote by z U) i, z V)N ^ the zeros of K v . It is known (cf. [TTJ) that the real 
part of each zero is negative and that the number N(y) is \u\ — 1/2 if v — 1/2 G Z 
and is the even number closest to \v\ — 1/2 if otherwise. 

Moreover, for /i ^ and c > 1, we set 

cos(tx^){I i1 (cx)K 1jL (x) - I^(x)K^(cx)} 
M ' cW " {K ll (x)} 2 + ^{I,(x)Y + 2nsm(ir^K fl (x)I ll (xy 

By the estimates for the modified Bessel functions, we see that L M;C (x) decay 
exponentially as x — > oo and that, 



(2.3) L„, c {x) 



l ° gC :{l + o(l)}, /i = 0, 



(log x) 
cos(7ru)(c^ - c M )a; 2M ri 

2^fWti) {l+o(l)l - " >0 - 



Theorem 2.1. Lei c > 1,1/ G R and w G C \ {0} wi£/i |arg(w)| < n and 
K v {w) £ 0. 

(1) J/i/ = ±1/2, we /tave 

i^(cw) _ e-^- 1 ^ 

(2) J/ |i/| < 3/2 and i/ ^ ±1/2, we have 

KJcw) e -^ w . n x... r°° we- {c ~ 1)x L 



K u (w) c\ u \ J x(x + w) 



X 
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(3) If v — 1/2 is an integer and v 7^ ±1/2, 



K v {cw) = e~^ w _ e(c _ 1)w K we^-V*^ K„(cz v>j ) 
(4) If u — 1/2 is not an integer and \u\ > 3/2, 



Jo x(x + w) 
Combining the results in Theorem 12.11 with formula 



CO 



e- xt q{t,a,b)dt = e~ (a - b)V ^, 



o 



where 

q(t, a, b) = _^i e -(a-6) 2 /2t a > b t > 

we can invert the Laplace transform (12. 2 j) and obtain an expression for the density 

» 

a, ft " 

We put c = a/6 > 1 and define the following functions. 



q(t,a,b), 



of r: «: 





#(0 f t 






iV(l/) 






= fl 






1 r 00 

" 6 Jo 



Zu,j 1 ( ^u,j , 

POO 



e (c ^[^(a^da; • q(t,a,b), 

L\ v \ >c (x)dx / £e~^~i?d£. 

y2iit 6 Ja-b 



Then we obtain the following from Theorem 12.11 

Theorem 2.2. TTie /irst hitting time has the density f^ which is given by 
the following. 

(1) If v = ±1/2, w;e have 

f^(t) = c-^q(t,a,b). 

(2) // < 3/2 and v ^ ±1/2, we have 
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(3) If v — 1/2 is an integer and v ^ ±1/2, 

/#(*) = c-^q(t,a,b) - c-^ a f(t) - c-^\t). 

(4) If v — 1/2 is not an integer and \u\ > 3/2, 

/£?(*) =c-^ H g(t,a,6) -c-^W -c-^ 2 (t) 

3. Asymptotic behavior of the densities 

In this section we study the asymptotic behavior of the density f^l of as 
t — > oo by using the expression given in Theorem I2.2L It is shown in [2] that 
/#(*) = c r 1 (logt)- 2 (l + o(l)) if z/ = and = c^- 1 "!-! (1 + (l)) if v ± . 

At first we note that, for every i/Gl, the functions a, 6), $i^' 2 (^) 
and ^it 1 ^) are 0(r 3 / 2 ) as t ->• oo. About ^ ,2 (t), since 



/■oo 

/ ^e- x ^ b d^ = 0(x- 2 ) as a; 4,0, 

«/ a— 6 



we obtain 



L\ v \ iC (x)—dx = oo 

from Q. This means t 3/2 ¥^ b ,2 (t) ^ oo, t ->• oo when < < 1/2. Hence 
we need to consider separately for the four cases, v = 0, < \u\ < 1/2, \u\ = 1/2 
and \u\ > 1/2. 

Proposition 3.1. Assume < b < a and put c = b/a. 

(1) If v — 0, we /lave 

(3.1) /^(t) = 21ogc-r 1 (logt)- 2 -(l + o(l)). 

(2) 7/0 < \v\ < 1/2, we /iaue 

(3-2) = ; 2Mr( ; D V "-(l + .(!))■ 

Remark 3.2. When > 1/2, we can prove that (J3~j2]) holds if v - 1/2 g Z. 
When I/-1/2 6Z and z/ 7^ ±1/2, we can prove that 

/2(i)^Mr [H - 1/2l - 3/2 '(i + o(i)) 

holds for some constant in a similar way to that in Section 4 of [4j , where we have 
studied the tail probability for rfl. We can give an expression for the constant 
C(u), but it is so complicated that we omit it. We believe that C(u) coincides 
with the constant on the right hand side of (13. 2p . In order to prove the above 
mentioned results when \v\ > 1/2, we need to use the result in [2] on the order 
of decay of the density. 
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In the rest of this article we directly deduce (13. ip and (13. 2p from the expression 
for f^(t) given in Theorem 12.21 

Proof of (pTTp . We have 

/$(*) = I 1 " ^e-( c -^L , c (x)rfx}g(t,a,6) + *<g' 3 (t), 
where c = a/b and (t) is given in the previous section. Setting 



S(t) 



'2Td 



a-b\ _i( n+ ^)2 



we have *i°l' 2 (t) = b~ l S(t). 

Fix arbitrary £ > 0. Then, by ( I2.3p . there exists 5 G (0, 1) such that 



log c 



(log x) 



< 



(logx) 



holds for every x G (0, 5). We also fix 7/ > 1 and assume that t satisfies 

Sa*f<s, „ < min { VMi ogt «, T i^ 5 <i 

V t log V t ^ 

We devide the integral defining S(t) into four parts, that is, we set 



Si(t) 

s 4 (t) 



1 
1 
1 



U + 



U + 



u + 



u + 



a-b\ _i( u+ £^)2 



Vt 



a — b\ _i[„i «=*i2 



e —2^—r du I L 0c (x)e- {c - 1)x e-^dx 



X y/tu 



Vi 



e'^Wdu 



Vi 



Vi 



L QtC (x)e~^ c ~ 1 ^ x e~^ :!L dx, 

(log t^M 

(logt)V^ 

L Ax)e-^e-^dx, 

v/Vt 



For an estimate for Si(t) we recall lo(-z) = 0(z 1 ^ 2 e z ) and iCo(z) = 0(z l l 2 e z ) 
as z — > oo. Then we see that there exists C\ > such that 



Hence we obtain 



ISiWI 



< 



00 / a-b 
u 



Vi 



du 



_0 xVtu 



which is 0(r 3/2 ). This implies t(logt) 2 S^t) 0, t -> oo. 
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By the choice of the constant 8, there exists C2 > such that 

C 2 



(3.3) 
Hence we get 



\Lq, c (x) 



< 



(log a:) 



< x < 8. 



< 



< 



<2nt 
C 3 



00 ' a-b 
u + 



y/i 



du 
dx 



(log t) 3 M ( lo g x ) 



y/tux 

e b dx 



t(logt) 3 7 ( iog t )3 /v ^a;(loga;) 2 ' 

From this we easily conclude that 5*2 (t) = 0(t _1 (logt)~ 3 ). 
For an estimate for Ss(t) we write 

^(logt) 3 /^ 



J2nt J \a-b)/ v 1 
Then, by ( I3.3j) . we get 



ue 



- uZ ' 2 du 



\/ / 1.; X 



V /Vt 

(togt) s M 
v/Vt 



Ln r (x)e b dx. 



1 y/tXU 

e 6 ctx 



(log 2; 



< 



C, 



t(logt) 2 ?7 



for some positive constants C4 and C5. 
We next set 



5 4 (t) 



1 



f 2id 



00 / a-b 
u 



ri/Vt 



Vi> du 



logc o _ {c _ 1)x Vt 



(logs) 5 



e b dx. 



Then we have \S±{t) — S±{t)\ ^ eS^{t) and 



taogVt) 2 ^^) 



log c 



~ u2 / 2 du 



'2lT J(a-b)/Vt 

By the choice of rj and i, we have for < y < r] 

log y/i , log y/i 



lo 



g y/t \ 2 -m 
— e b 



\log(Vi/y)J 



dy. 



< 



]og(y/i/y) log y/i- log T] 
Hence the dominated convergence theorem implies 

li mt (log^ l(t ) = ^-^/ 

t^oo 2 yJllX Jo 



< 2. 



e 2 6 du. 



Combining the above mentioned estimates and letting rj — > 00, we obtain 

b logc 



lim sup 

t— >oo 



t(log \/t) 2 S(t) 



, 6 logc 
^ — - — £ 



and the desired result. 



□ 
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Remark 3.3. By using the formula (see [HI p. 80]) 

K (z) = - \og(z/2)I (z) + J2 ~T^^ m + *)> 

where ^(m + 1) = Y^k=\ k~ l — 7 for the Euler constant 7, we have asymptotic 
expansion for L 0>c (x) as x j. and can give the asymptotic expansion for ff$(t) 
in the form 

(0) 21ogc f «!_ Q2 . . 

Ja - fel ; t(logt) 2 \ log* (logt) 2 

After some computations, we can show ot\ = 2(7 — log 2 + 2 log b). In general 
the explicit expressions for the constants 0:2,03, ... are complicated and we omit 
them. For a related result, see fTUj. 

Proof of We start from 

/£?(*) = c-^(*,a,6){c-M - h H>c (x)e- {c - 1)x dx\ + C -^ b \t). 

Then, by a simple change of variables in the defining integral for $[j 2 (t), we 
obtain 



6V27rt 1 +l"l 



00 ^00 



= ; 7=^ / / l[(a-6)M,oo)(«) 



</0 



(y/v^)2|H 

By (12. 3p the function L\ v \ >c (x) / 'x 2 \ u \ is bounded near and, on (0, 00) since L\„\ )C (x) 
decays exponentially at 00. Therefore, by the dominated convergence theorem, 
we obtain after some manupulations 



(3.4) lim t 1+ M*M' 2 (t) = -^=a,r(l + 2\u\)0 v h-^- 1,2 v( 

V27T V 



1 - 2IH 



where the constant C v is given by 



L Kc (x) C H_ C -H 
Gj, = hm — — = cos(z/7r) 



xio xW K 7 2 2 IH-ir(|z/|)r(|z/| + 1)' 



Finally, by simplifying the right hand side of ( 13 ,4p via the functional equality 
and duplication formula for the Gamma function, we arrive at the desired result. 
□ 
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